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The density matrix expansion is used to derive a local energy density functional 
for finite range interactions with a realistic meson exchange structure. Exchange 
contributions are treated in a local momentum approximation. A generalized Slater 
approximation is used for the density matrix where an effective local Fermi mo- 
mentum is chosen such that the next to leading order off-diagonal term is canceled. 
ffartree-Fock equations are derived incorporating the momentum structure of the 
underlying finite range interaction. For applications a density dependent effective 
interaction is determined from a G-matrix which is renormalized such that the 
saturation properties of symmetric nuclear matter are reproduced. Intending appli- 
cations to systems far off stability special attention is paid to the low density regime 
and asymmetric nuclear matter. Results are compared to predictions obtained from 
Skyrme interactions. The ground state properties of stable nuclei are well repro- 
duced without further adjustments of parameters. The potential of the approach is 
further exemplified in calculations for A=fOO. . .f40 tin isotopes. Rather extended 
neutron skins are found beyond ^"^"Sn corresponding to solid layers of neutron mat- 
ter surrounding a core of normal composition. 
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I. INTRODUCTION 

Studies of nuclei far off stability have given a new impact to nuclear structure physics. Exotic 
nuclei allow for the first time to test nuclear theory at extreme isospin and, close to the driplines, 
also at low density. In the past nuclear models have been derived almost exclusively from well 
bound and charge symmetric stable nuclei. Hence, structure investigations far off stability are 
providing important extensions into hitherto unexplored regions of the nuclear landscape. It is 
apparent that such studies involve far reaching extrapolations of nuclear models much beyond the 
range where they were initially developed. The weaker binding of strongly asymmetric systems 
enhances interactions between bound and unbound configurations and the regime of open quantum 
systems might be approached. In extreme dripline nuclei as e.g. ®B Q] and ^^Li |^ evidence for a 
strong competition of interactions between individual nucleons and mean-field dynamics is found. 

At present nuclear structure theory is making first exploratory steps into the new regions far 
off /9-stability. Experience shows that seemingly well established models describing stable nuclei 
equally well are leading to strong deviations when applied to exotic systems. A particular striking 
feature are the variations in predictions of binding energies and other ground state properties 
at large charge asymmetry. The results indicate that the isovector properties of interactions are 
not well enough understood in order to reproduce or predict the location of driplines. A widely 
applied method in non-relativistic structure calculations is the Skyrme approach Although 
very successful in stable nuclei, it is an open question to what extent the restricted set of operators 
used in the Skyrme energy functional is sufficiently accounting for the conditions in the new mass 
regions. The transparent structure and the good applicability is a clear advantage of the Skyrme 
model. It is also fiexible enough for adjusting the already incorporated parameters or introducing 
additional ones in order to describe new data. However, a purely empirical approach might be 
hampered by increasing uncertainties when approaching the limits of stability. 

The intention of the present work is to indicate an alternative approach which is closer to a more 
microscopic description but still retains to a large extent the advantages of using a local energy 
functional as in Skyrme theory. A more fundamental description requires to use microscopically 
derived in-medium interactions as e.g. the Brueckner G-matrix accounting for medium effects 
on the level of ladder diagrams Theoretically, such a full scale many-body calculation is 

not feasible as a standard approach to finite nuclei, except for selected cases A successful 

program is to use instead a G-matrix from infinite nuclear matter and to apply it in the local density 
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DTOximation (LDA) to finite nuclei [Q. The LDA has been used succesfully in non-relativistic 
and also covariant fieldtheoretical models [jllJ-Q. However, the finite range character of 
a realistic G-matrix is well known to introduce complications on the numerical level by leading 
to a system of integro-differential equations |]l9| . A way around this problem was discussed some 
time ago by Negele Q but rarely applied later. The numerical complications arise because of the 
non-locality of exchange contributions to the HF-field. For a finite range interaction they involve 
the full one-body density matrix rather than the local densities needed for the direct parts and 
contact interactions of Skyrme type. Negele introduced an operator-valued expansion of the non- 
local parts of the density matrix. The contraction with the HF-ground state configuration leads to 
a systematic expansion in terms of the local one-body density. The non-localities are contained in 
functions where the lowest order contribution coincides with the well known Slater approximation 
1^,^. This density matrix expansion (DME) allows to integrate out the non-localities of the 
exchange terms. 

The approach presented below follows in spirit closely the original work of Ref. §,§. In section 
the basic steps of the DME are summarized. Different to Negele we choose an effective local "Fermi" 
momentum q^{r) such that the next to leading order term of the DME series is canceled exactly. 
q|n is determined by the kinetic energy density and, in addition, the Laplacian of the one-body 
density. For an infinite system the conventional expression for the Fermi momentum is recovered. 
The DME results are then used to define a local energy density functional. Since it originates 
from an effective interaction, assumed to be given independently e.g. by Brueckner calculation, the 
parameters of the model are fixed. By variation of the energy functional HF-equations are derived. 
With our choice for qp the kinetic energy part of the single particle Schroedinger equation receives 
additional contributions. Thus, we obtain consistently effective masses m*„(r) for protons and 
neutrons, respectively, which are completely determined by the momentum and isospin structure 
of the original interaction. An important extension over the Skyrme approach is that the self- 
consistent local momentum qp contributes in arbitrarily high order to the exchange terms which 
is especially important at high and low densities. 



Applications to infinite symmetric and asymmetric nuclear matter are discussed in section III 



We start by choosing a somewhat simplified effective interaction. The three Yukawa parameter- 
ization of a G-matrix calculated initially for ^^O by Toki et al. 22| is used. This so-called 



M3Y interaction has been widely and successfully applied in nuclear reaction [g3|,^ and struc- 
ture calculations |p^ , p5| , p6| . The former studies, especially for nuclear reactions, indicate that this 
interaction has a realistic momentum structure including a long range one pion-exchange tail. It 
is especially trustable at low nuclear densities around 1/3 of the saturation density. A density 
dependence is lacking and it is known that the original M3Y parameterization does not lead to 
saturation of nuclear matter. As other authors before |24|, we restore the density dependence 
by fitting the properties of symmetric nuclear matter at saturation including the binding energy, 
equilibrium density, compressibility and, in addition, also the symmetry energy. This allows to 
fix the density dependence of the isoscalar and isovector vertices in the spin-scalar particle-hole 
interaction channels. We thus have determined a renormalized density dependent three Yukawa 
(D3Y) interaction which we expect to resemble closely the momentum structure of a G-matrix but 
with an effective strength fitted to nuclear matter properties . The in-medium particle-particle 
interaction is obtained by the appropriate transformation and used in pairing calculations. At low 
densities the strength of the free singlet-even NN T-matrix is approached. Results for symmetric 
and asymmetric nuclear matter are compared to calculations with empirical Skyrme interactions. 
Close to saturation a good agreement is found but deviations occur away from that point. 



In section IV the model is applied to finite nuclei. The parameters are kept fixed as ob- 



tained from the nuclear matter fits. Binding energies, charge densities and radii of stable nuclei 
^i6Q^4o,48Qg^^902j- 208p|-,-j g^j.^ ^^^^i described. The good agreement confirms the approach and 
makes applications to unstable nuclei meaningful. HF calculations for Sn isotopes lead to a good 
description of the binding energies for A>120 but a slightly too weak binding is found in the 
lighter nuclides. The calculations predict rather thick neutron skins for A>130. The paper closes 
in section with a summary and conclusions. 
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II. DENSITY-MATRIX EXPANSION AND HARTREE-FOCK THEORY 



The total Hartree-Fock energy of a nucleus with A nucleons and the single particle states |fc) is 
given by the kinetic energy 



k<A 



(r)- / d\— {Tp(r)+T„(r)}, 



(1) 



where q — n,p denotes protons and neutrons and Tg(r) the kinetic energy densities. The potential 
energy of a two-body interaction V 



{V) 



Y,{k^k2\V\k^k2) - (fcifcal^lfefci) 

.felfe2 



(2) 



is expressed in terms of interactions for like {q = q') and unlike (q ^ q') particles as defined in 
App. 

{V) - iV) + {V^) = E {K') + iK,')} (3) 

qq'—p^7i 

= E/ ^'''^1 / '^^'^2 {p,(ri)p,,(r2)F/,,(ri2)+p,(ri,r2W(ri,r2)F/,,(ri2)} 

qq' 

where direct and exchange contributions are indicated separately by the index d resp. e. Here, 
ri2 = ri — r2 is the relative coordinate and the density matrices are defined as 



Pqiri,r2) = ^0fc(ri,cr,q)0fe(r2,cr,q) 

A; (T 

where (r2 , cr, g) denotes single particle wave functions. 



(4) 



A. Density Matrix Expansion 

The exact treatment of the exchange term in the Hartree-Fock functional would, after variation, 
lead to a coupled system of integro-differential equations. A numerically simplified and theoretically 
more transparent approach is obtained using the density-matrix expansion (DME) which was 
originally invented by Negele and Vautherin ||^. In the DME approach one first transforms to 
center-of-mass and relative coordinates by introducing r = ^^^j^ and s = r2 — ri. Then the density 
matrix is formally expanded into a Taylor series in s. A particular simple result is obtained if the 
time-reversed orbitals in the nucleus are filled pairwise g]. Therefore, the DME approach is most 
suitable for even-even nuclei. As shown in [^, averaging over spin directions and resummation 
leads to a series expansion where the first two terms are 

p(r+^,r-^) = psL{sqF)p{r) (5) 



2' 2 



35 

—3-j3[sqF) 



2sq-y 



iv2p(r)-T(r)-K^4p(r) 



and in the leading order term the well known Slater approximation |^,^ is recovered. Here, the 
Slater density psLisqp) = ■^^jiisqp) is given by the spherical Bessel function of the order 

n = I. 

The quality of the expansion depends on the choice of q|, (r) which should obviously correspond 
to the average relative momentum between the two interacting particles. Different to Ref. we 
choose 

5 r(r)-lvV(r) 
'^-(^) = 3 ■ 
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which is seen to cancel exactly the next to leading order term in Eq.(|5|) and has the property 
Jd^rp{r)ql{r) = 5/3 J d'^r T{r). 

With this choice the DME reduces to the well known Slater approximation /3(r+|,r — |) — 
PsLisqp) p(r) of the exchange term but now with a modified Fermi momentum that accounts 
for surface corrections to the Fermi momentum of the Thomas Fermi model up to the second 
order. The advantage of this method is that the description of the exchange potential is simplified 
and can be treated in a way similar to the nuclear matter calculation. It is easily seen that our 
choice of the Fermi momentum reduces in infinite nuclear matter to the conventional expression 
kp = (3/27rV)'/^. 

Once the relative and the center-of-mass coordinates are separated and the integration over 
the relative coordinate s can be immediately performed analytically. The potential energy of the 
exchange part finally simplifies to: 

{Vn=Yl J d\pq{r)p,,{r)V^^,{r) (7) 

qq' 

In this equation the strength of the interactions are defined as 

Kq,iqF,,qF^,;r) = J dhV^q,is)psL{sqF^{r))psLisqF^,ir)) (8) 

In this way the exchange term of a finite range interaction is mapped to an effective contact 
interaction with a local strength depending on the effective local Fermi momentum, Eq. (^) , at the 
place of the interaction 

V^\, is, r) ^ SHs) V^,, [qF, (r), g^,, (r)). (9) 
We are now at the position to reformulate the theory in terms of a local energy density functional 



with the total energy 



H{v) = Ht{y) + Hv^{v) + Hv^r) 
E = j d^rH{v) = r + i(t^). 

utq'{v)^ j dV p,,{v')v^%,{vy) 

Ulq'{v)^Pq'{v)V^,.{q\,q\,-.v) 

B. The Hartree— Fock Equations 



(10) 

(11) 

(12) 
(13) 

(14) 

(15) 
(16) 



Comparing this to the equations derived before one finds: 

i?T(r) = |l [r,(r)+r„(r)] 

qq' 

Hv4v)^\Y.Pq{v)Ul^,{v) 

qq' 

Single particle potentials U{r,p) have been introduced as e.g. 



In our case, the total energy depends in a complex, nonlinear way on the nuclear density p, 
the kinetic energy density r and V^p. By using the DME we are able to simplify the exchange 
term. Our choice of the Fermi momentum qpip, t, V^p), Eq.(||), leads to an additional functional 
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dependence on the kinetic energy and the Laplacian of the density and therefore takes into account 
dynamical and surface corrections. For the variation of the HF energy functional, Eg. (|lo|) , we 
express p, r and qp by 0^ and its derivatives: 

p{r)=Y.\Mr)f (17) 

k 

T(r)=5]|V0fc(r)|^ (18) 

k 



The HF variational equation 

" '^\Mr)\^^k] =0 (20) 



J d^rS (^H{r)-Y^ 



leads after partial integration to 

.e.i:/a..,.(|^-v^-v^j^-*«) (21) 

For the further calculations we assume that the nucleon-nucleon interaction has an intrinsic 
density dependence resulting from the medium dependence of the meson-nucleon vertex that can 
be parameterized by vertex corrections q-y{p). A model will be specified in the next section. 

Then, the partial derivative in Eq.(|2l|) with respect to (t)k is written as 

dH _ dH dp dH dg^jp) dp dH dql 

d^l dpd^k dg-^ip) dp d4>l dqld^k ^ ' 

The first term is the ordinary HF potential, the other terms are rearrangement terms resulting 
from the intrinsic density dependence of the interaction. Analogously the other partial derivatives 
are expressed as 



dH _ dH dql 



dH dH dr dH dql 

^7 (23) 

k 

(24) 



9V2(/.* dqj, dy^(t)i 



Here, the second term in Eq. d23h g ives rise to an effective mass. Using the above definitions and 
plugging in everything into Eq.(|21|) leads to the Schroedinger-equation 



2m* (r) 

The Hartree-Fock potential is given by 



V + U^' (r) + U^{r) 0,,(r) = e,,0fe,(r). (25) 



'^r(r)=E[<'W + f^.V(r)] (26) 

and the effective masses are 



= r::- ^ • (27) 



1 _i_ ^ 



The effective mass results from the dependence of the Fermi momentum on the kinetic energy 
density r in the exchange term. It is a direct consequence of our choice of the Fermi momentum 
and the finite range nature of the interaction. 
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The variation of the density dependent vertex functions g-yi^p) and of the Fermi momentum in 
the Fock energy density leads to an additional rearrangement potential 



' ^ ' dg,{p) dp 



(28) 



1 2 



dql dql 



24 



"^qq I 2 W 



dql dql 



which is of importance for the mean-field or single particle potential in the HF equations. It results 
from the density dependence of the NN interaction and is known to describe (static) polarization 
of the background medium |l8J2q,M. 



III. NUCLEAR MATTER 
A. The Nucleon— Nucleon Interaction 

A well established approach to derive in-medium interactions from the NN-interaction in free 
space is Brueckner theory. However, the ladder approximation, even if refined by going beyond the 
quasiparticle approximation as e.g. in Ref. |]30| , is generally found to miss the empirical saturation 
properties of nuclear matter at least when including two-body correlations only. For the purpose 
of the present work we therefore use a semi-phenomenological approach. 

We start from the M3Y G-matrix [|l|,|2) obtained originally for ^^O. The M3Y G-matrix is 
parameterized by a superposition of three Yukawa functions where the ranges were chosen to 
represent the long-range tail of the one-pion-exchange potential (OPEP) and medium and short- 
range parts. The shorter ranges correspond to boson masses m = 490 MeV and m = 780 MeV 
and thus mimic a, lj and p meson exchange interactions. Here, we use the parameter set obtained 
from the Paris NN-potential j2^ which is given in Tab. |. 

Because of the lack of a density dependence the M3Y interaction can be considered as realistic 



only for a small range of nuclear densities |31| being centered around one third of the saturation 
density. In connection with optical potential and nuclear reaction studies various authors, e.g. 
Kobos et al. and more recently Khoa and von Oertzen |^^, have noted that the "bare" 
M3Y interaction leads at high densities to the collapse of nuclear matter but by supplementing a 
phenomenological density dependence the correct binding properties could be restored. In spirit 
we apply the same approach by introducing density dependent vertices fitted to the saturation 
properties of infinite symmetric nuclear matter. 

In the former works only the medium dependence of the isoscalar-spin scalar channel was con- 
sidered [^4|j32[ . Here, we choose a more general description and include a density dependence also 
in the isovector channel. We introduce density dependent isoscalar and isovector vertex functions 
fo,T(p), respectively, multiplying the corresponding "bare" M3Y vertices as a whole but leaving the 
intrinsic momentum structure untouched. Hence, we assume that the many-body and the momen- 
tum structures are separable in which we follow the observations made in determining effective 
in-medium meson-nucleon coupling constants from Brueckner calculations p^- p^ , |33[ | . 

Only central interactions are considered. Aiming at applications to asymmetric systems it is 
suitable to rewrite the central singlet and triplet even and odd components [Vse, Vte, Vso, Vto) 
in terms of interactions for like (proton-proton, neutron-neutron) and unlike (neutron-proton) 
particles. We then consider the interaction in the "particle-hole" channel which eventually deter- 
mines the Hartree-Fock mean-field. Multiplying the spin and isospin exchange operators onto the 
interaction effective direct (d) and exchange (e) interactions are obtained Q|.With the relations 
given in App. ^ the isoscalar and isovector interactions are then obtained as 

= I + + + Vnl') (29) 

yd/e ^l^yd/e _yd/e _yd/e ,yd/e^ /gQN 

*r PP pn np ' ^ nn / J 

and medium dependent vertex renormalizations fo,T(/o) are introduced 
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y*,'J\p) = foApKi' (31) 

taken to be the same in the direct and the exchange channel. Back transformation to the sin- 
glet/triplet representation (see App. ^ is used to define the corresponding in-medium particle- 
particle interaction. By this procedure we retain the basic relations of particle-hole and particle- 
particle interactions. 

The question arises how to choose the functional form of the vertices. As a guidance we refer 
to Brueckner theory and express the G-matrix R in terms of the free NN scattering matrix K and 
subtract the Pauli-blocked contributions from inside the Fermi sphere p^ , p5[ |: 

Rik, k') = K{k, k')-P [ d^qK{k, q)GF{q)R{q, k') (32) 



where Gp denotes the two-body propagator inside the Fermi sphere. Dependencies on the two 
particle center-of-mass momentum and energy have been omitted. For the coupled partial wave 
channels Eq.(|3|) leads to a linear system of equations. Without attempting to solve Eq.(|3^) exactly 
an estimate of the density dependence is obtained by choosing k=ki? and, in addition, also q=kir 
in the interactions beneath the integral. This amounts to the reasonable assumption that the 
momentum dependence of K and R over the Fermi sphere is much weaker than the variation of 
the (singular) quantity Gp- The principal value integral over Gp can be performed analytically 
and one finds that the medium corrections are given in leading order by the density of states at 
the Fermi surface, < Gp >~Ni?=mki?/(7r?i)^. (< Gp > depends in fact also on the two-particle 
center-of-mass momentum and energy which we assume to be fixed to appropriate average values). 
Thus, as an approximation to Eq.(p2|) we find for the half-off shell G-matrix 



Next, we consider Eq.(32) for the case that a solution Rq at a Fermi momentum k^ is known. 
In Eq.(p3|) this corresponds to insert in the denominator the subtracted Green function < Gp > 
— < G^f>. Since G^ is a constant and K is weakly varying we may rewrite Eq.(|3^) as 

R{kp, k) ^ ,.o^p rt.^ . wr^ s{ky)R^{kp, k) (34) 
1 -I- s(k^f)Ro[kp, kp) < Gp > 

and a scaling factor 

'(^^^ ^ l-i?o(fc°,fc°)<G^> ^^^^ 

is found. 

These considerations suggest to use for the vertex functions (7 — 0, t) the ansatz 

N-, 

Mp)^s,\l + Y,aZz^^\ (36) 

n=l 

where z=p/po and the base exponent is fixed to P=l/3. The coefficients a^ are to be determined by 
the fit to nuclear matter properties. As seen below nuclear matter allows a first (Nt- = 1) and third 
(No =3) order approximation for f,- and fo, respectively. The overall scaling factors s^ renormalize 
the interaction strength asymptotically for p — > such that with an appropriate readjustment of the 
a,n coefficients the nuclear matter saturation properties are unaffected. From the above discussion 
it is seen that the scaling factors relate the interaction vertices to the initial renormalization point 
k^. In the following, s-y and a^ are treated as empirical parameters. 
The ground state energy of nuclear matter is 



E E I? + J E E E a'k'l V |ak, a'k') (37) 



. _ ,^2^,2 

E , , 

2m 2 

a k<kF OLOi' k<kF k' <k' 
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where V includes antisymmetrization as discussed above and a represents the spin an isospin 
quantum numbers a and q. In nuclear matter the single particle wave functions \ak) are plane waves 
and the Fermi momenta kp^ of protons and neutrons, respectively, is related to the corresponding 
densities pq by kp — (Stt^p^)^/^. The total binding energy per particle in nuclear matter is 
obtained as: 



E{p,Pt) _ 3 TpPp + TnPn 

A ~ 5 p 

qq' 



(38) 



Here, p 



I + pp and Pt = Pn — Pp are the isoscalar and isovector density, respectively, and 

iq — denotes the Fermi energy of the protons and neutrons. The effective density dependent 

interactions appearing in the above expression are given in nuclear matter as the volume integral 
of V*'^{s,p) for the Hartree contributions, e.g. 



qq' 



(39) 



while the non-local Fock contributions are averaged over the Fermi spheres, e.g. 

Kq'ip;f^F,,kF,) ^ / <fspsLiskFjpsLiskF,)V*5{s,p) 



(40) 



For Yukawa form factors the exchange integrals can be evaluated in closed form which, however, 
is omitted here. In symmetric nuclear matter, Eq.(B8h simplifies to p0|,p4| 



(41) 



|(P, = 0) = + i p [v,'{p) + Vo'ip, kp) 



^T+^pfoip) 



where kp = (3/27r2p)i/3 and r = For given So,T and exponent P the remaining parameters of 

fo are determined by imposing the conditions for the nuclear matter binding energy, the saturation 
density and the compressibility at the equilibrium density po, 



(p, pr = 0) 



E 
A 



2 diE/A) 



9p 



dp 



2 d^E/A) 



-ay 
= 



dp2 



(43) 
(44) 
(45) 



and the isovector vertex renormalization is obtained from the symmetry energy e-r{p), defined 

by 



E{p,pr) Eip,0) 



A 



A 



St{p) 



£4(p) 



(46) 



At the saturation density p ^ po the first term in Eq.(|4q) gives the binding energy per particle of 
symmetric nuclear matter, —ay, and Et is the symmetry energy coefficient 



2 d^E{pQ,pr) 
as = Po 



dpi 



(47) 



Px=0 



which is known to range between 28 MeV and 32 MeV |20 3£ 
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For a finite range interaction the above procedure for fixing tfie isovector properties of asymmetric 
nuclear matter involves higher order contributions from the intrinsic momentum structure of the 
interaction. They appear because the asymmetry in the proton and neutron Fermi momenta 
induces an effective isospin symmetry breaking component via the exchange integrals extending 
over the two Fermi spheres. Numerically, the mixing of the isoscalar and the isovector channel is 
noticed immediately from the separate averaging of the exchange contributions over the proton and 
neutron Slater densities psLikpp) and psLikpn), Eg. (pl|) . In the direct part or for a momentum 
independent contact interaction, respectively, these contributions are absent. In those cases the 
potential energy can always be written in the form p'^Vo+p^Vr, which separates exactly the isoscalar 
and the isovector part of the interaction (except for the case that Vq would depend intrinsically 
also on p-r). Skyrme interactions account to some extent for the isovector exchange effect by the 
momentum dependent components. At extreme high asymmetries, however, differences to a finite 
range interaction might become noticeable. 



TABLE I. Parameters of the unrenormalized M3Y interaction in the odd-even (first block), the 
spin-isospin (second block), the p-n representation as defined (third block) and the effective spin-isospin 
representation (fourth block) as defined in App. 



ro [fm] 


0.250 


0.400 


1.414 


Volume integral 


m=^ [MeV] 


780 


490 


140 


[MeV fm^] 


VsE 


11466.00 


-3556.00 


-10.46 


-980.28 


Vte 


13967.00 


-4594.00 


-10.46 


-1324.02 


Vso 


-1418.00 


950.00 


31.39 


1600.77 


Vto 


11345.00 


-1900.00 


3.49 


823.43 


Vo 


11062.00 


-2538.00 


0.00 


131.17 


14 


939.00 


-36.00 


0.00 


155.40 




314.00 


224.00 


0.00 


241.33 


Var 


-969.00 


450.00 


3.49 


295.53 


yd 


11375.00 


-2314.00 


0.00 


372.51 


yd 


10748.00 


-2761.00 


0.00 


-110.16 


ypp 


-5642.00 


536.00 


-5.23 


-862.65 


^ pn 


2594.00 


-1574.00 


-10.46 


-1127.93 




11062.00 


-2538.00 


0.00 


131.17 




314.00 


224.00 


0.00 


241.33 




-1524.00 


-519.00 


-7.85 


-995.29 




-4118.00 


1055.00 


2.62 


132.64 
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B. Results for Infinite Nuclear Matter 



The isoscalar vertex function io{p) was adjusted to the empirical vahies of the binding energy 
(ay ~ 16 MeV), the saturation density {po ~0.159 fm"'^) and the compressibiHty K^o- Empirically, 
Koo is not well known (values range from 200 MeV to 350 MeV), but recent studies, e.g. Q, 
converge to a value of about K(x>=230 MeV. For the symmetry energy coefficient serving to fix ir 
as = £^(po=0.159 fm-3)=32 MeV is used. 

In Tab. || two sets of vertex renormalization parameters, Eq.(36) are displayed. The first set 
was obtained with Sq = 1, i.e. leaving the free space interaction unchanged, whereas the second 
set with So = 5 rcnormalizcs the asymptotic behavior of the interaction as well. This value of Sq 
was determined from the binding energies of finite nuclei as will be discussed in the next section. 



TABLE II. Parameterization of fo(p) and irip) 







al 






po [im-'-"] 


K [MeV] 


{Eg /A) [MeV] 


So 


= 1 


1.365 


-1.839 


0.379 


0.159 


230 


-16.00 


So 


= 5 


-2.127 


2.032 


-0.724 


0.159 


230 


-16.00 






al 










er{po)[MeV] 


St 


= 1 


.1688 










32.00 
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1 1 1 1 1 1 1 1 1 1 1 1 1 1 

0,0 0,1 0,2 0,3 

nuclear density p [fm"^] 

FIG. 1. Density dependent isoscalar fo (so — 1 dashed line, sq = 5 solid line) and isovector fr (dotted 
line) vertex functions of the D3Y interaction. 

Fig. 1^ shows the density dependence of the vertex correction functions fg and f^-. In the isovector 
channel only a rather small deviation from unity is found, and f,- increases only slightly for higher 
densities. This reflects the fact that the original interaction already has a symmetry coefficient 
close to 32 MeV. In the isoscalar channel one realizes an increase of the interaction at low densities 
and a suppression of the interaction at higher densities. The low density behavior obviously scales 
with Sq. In Fig. || the equations of state (EoS) for the D3Y and the Skyrme SLy4 interaction 
are compared. For So=l the D3Y and the SLy4 results are in very close agreement over the full 
range of densities. According to the behavior of fo, the binding energy increases at densities lower 
and decreases at densities higher than the saturation density for Sq >1 while retaining its value at 
the equilibrium density. 
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FIG. 2. Nuclear matter equation of state calculated with the density dependent D3Y interaction for 
the parameterizations so = 1 (dashed line) and sq = 5 (solid line) and the Skyrme SLy4 interaction (dotted 
line). 



The isovector properties of the D3Y interactions are illustrated in Fig. where the EoS for 
asymmetric nuclear matter is displayed for different asymmetry ratios Z j A — Ppj p. Stable nuclei 
typically cover the range ZjA = 0.5... 0.4 (e.g. '"^Ca, ^"^Pb). In exotic nuclei the new regions 
with ZjA — 0.4 . . . 0.3 become accessible and far of stability ZjA = 0.2 may be reached (e.g -'^°He). 
It is remarkable that neutron matter is still bound but this occurs at extremely low densities. The 
reason for this behavior is the strong attraction of the D3Y with Sq = 5 at very low densities and, 
as such, is not very conclusive. In stable nuclei these density regions will rarely contribute to the 
total binding. However, a different situation may be expected far off stability, especially in halo 
nuclei, where the binding sensitively depends on the interactions at low density. Comparing the 
D3Y to the SLy4 interaction one sees that both interactions are in agreement around the saturation 
density for asymmetry ratios of ZjA = 0.5 . . . 0.4. The differences showing up at low densities and 
large asymmetries are mainly due to the stronger attraction of the D3Y isoscalar interaction with 
so=5. 
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FIG. 3. Equation of state for asymmetric nuclear matter calculated with the density dependent D3Y 
interaction for the parameterizations sq = 5 (solid line) and the Skyrme SLy4 (dashed line). Shown is the 
EoS for neutron matter and for the asymmetry ratios ZjA = 0.2, 0.4, 0.5. 

In Fig. ^ the EoS for different asymmetry ratios and the two sets of parameters are examined. 
The equilibrium density po(0' the energy per nucleon _Bb(^)/A, the symmetrie energy er(C) ^tnd 
the compressibility Koo(0 for i^^Z/A are shown. A higher value of So is seen to lead to a lower 
equilibrium density for strong asymmetries. In addition the binding energy at the equilibrium 
density is higher due to the stronger attraction of this parameterization of the interaction at low 
densities. The lower left part shows the symmetry energy of the EoS in dependence of the proton 
excess. One realizes a decrease of the symmetry energy for strong asymmetries and a higher value 
of So corresponding to the shift of the saturation density to lower densities. The compressibility 
approaches zero at high asymmetries with a remarkable flat slope. 
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FIG. 4. Equilibrium density po of the EoS (upper left), binding energy per nucleoli Eb/A at po (upper 
right), symmetry coefficient e,- at po (lower left) and the compressibility Koo(po) (lower right ), all plotted 
as a function of the neutron excess or asymmetry ratio ^ = Z/A. Results are shown for the D3Y interaction 
with So = 1 and so — 5. 



The properties of the D3Y interaction in the particle-particle channel were investigated by consid- 
ering like-particle pairing in symmetric nuclear matter. Pairing, being mediated by the singlet-even 
part VsE of the interaction, was examined by Hartree-Fock-Bogohubov (HFB) calculations. The 
in-medium SE particle-particle interaction was obtained by back transformation from the renor- 
malized particle-hole interactions using the representation in terms of the exchange interactions 
(see App. 0), leading to V^^ip) = 2E^=o,r [V^'ip) + V;%p)] ■ 

As in refs. [p7 38 1 the following set of equations was solved simultaneously: 



pikp) 
w^(fc, kp) 

Aik,kF) 



_fk_ 
(2^ 



w^(fc, kp) 

e{k, kp) — X 



2 I ^ ^{e{k,kp)-X)^ + A^k,kp) 



= 1 — u^{k, kp) 



2l2 



d^k' 



2m 'J (27r)3 
d^k 



{kk'\V*\kk' - k'k}v^{k',kp) 



{k - k\V^p\k' - k')uik\ kp)v{k', kp) 



(48) 
(49) 

(50) 
(51) 
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FIG. 5. Scattering length aFg^ikp) for relative s-wave states. Results are shown for the bare M3Y 
(dotted) and the D3Y with the parameterizations so = 1 (dashed line) and so = 5 (solid line). 



Fig. H shows the calculated in-medium scattering length for like particles (pp/nn) alg'^{kp) at 
a given Fermi momentum kp for different parameterizations of the interaction. Because of the 
density dependence in the isovector and in the isoscalar channel the strength of the interaction 
increases at low to intermediate densities. For the D3Y with Sq = 5 the pairing interaction is 
enhanced by a factor of about 8 at low densities over the values of the original M3Y because of 
the asymptotic scaling of the isoscalar interaction. This corresponds to about 90% of the s-wave 
strength of the free NN T-matrix with a scattering lenght of aFg^ = — 17.1fm and leads to higher 
values of the gap in nuclear matter. In Fig. ^ A(fc, kp) at k—kp is shown. For the bare M3Y 
the position of the maximum gap lies at about kp ~0.8 fm~^. About the same result is obtained 
in HFB calculations with the Gogny force |3^,Q and for several Skyrme interactions |3|]. The 
position of the maximum gap in the D3Y-calculations is shifted to slightly lower densities with 
kp ~=0.75 fm~^ resp. kp ~0.63 fm~^ for So=l resp. So=5. One finds that the height of the 
maximum increases by lowering the position of the maximum. Our results of A{kp, kp) ~ 2 MeV 
resp. A{kp, kp) ~ 4 MeV for the M3Y and the D3Y with So = 1 are in agreement with the Gogny 
force that yields a value of A{kp, kp) ~ 3 MeV |37 3^. The calculated maximum gap for the D3Y 
with So = 5 has a much higher value of about 9.5 MeV. The reason for this is the much stronger 
pairing interaction, especially at lower densities, as seen in Fig. ^. Comparing the value of the 
gap at the saturation density kp = 1.33/m~^ with other calculations we find a value of about 0.2 
MeV for the M3Y and a nearly vanishing value for the D3Y interaction. For the Gogny force one 
finds a value of about 0.6 MeV but most other calculations favor a small to vanishing gap 
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FIG. 6. Gap at the Fermi momentum in nuclear matter. Results for the D3Y interaction with sq = 1 
(dotted line) and sq = 5 (solid line) and the Gogny Dl (dashed line) are shown. 



IV. HARTREE-FOCK CALCULATIONS FOR FINITE NUCLEI 



The density dependent D3Y interaction derived in the preceding section for nuclear matter is 
applied in HF-calculations for finite nuclei. Stable nuclei close to the N~Z line are chosen first. 
In these calculations mainly the isoscalar properties of the interaction are investigated. A more 
general test, including also isovector dynamics, is obtained from calculations for even-even Sn- 
isotopes covering the range from the proton dripline (^°"Sn) to the A=140 mass region. The 
calculations for stable and unstable nulcei are performed without readjustment of parameters from 
their values derived for nuclear matter. 

The general structure of the HF-equ ation s for a density dependent NN-interaction in the DME 



approach was already derived in Sect. II B . In the practical applications the DME is used only 
for the exchange contributions to the HP potentials. Especially at higher density and thus Fermi 
momentum exchange is dominated by the short range pieces of the interaction for which the local 
momentum approximation underlying the DME is most appropriate. Different from the original 
work of Negele and Vautherin the direct, i.e. Hartree potentials are treated in full finite range. 
The D3Y density dependence was introduced as a local renormalization of the vertices which 



requires to replace the interaction by |18 



V;\v^,V2;p)^9,{v^)V^{v,~V2)g^{v2) (52) 

with = /o,r and where ^^(r) = g^{p{Y)) are the in-medium isoscalar and isovector vertex 
renormalizations, respectively. V^{yi — Y2) denotes the unrenormalized "bare" interaction in these 
channels. This ensures a reliable treatment of the density dependence also in the surface region 
where the vertices vary rapidly (see Fig. The direct potentials are obtained as in Eq.(^5|) 
by folding the proton and neutron densities, respectively, with the corresponding - now density 
dependent - direct interactions. To accomplish this we express the interactions V]j^/(ri,r2) by 
V^*''(ri, r2; p). With the DME the exchange or Fock potentials are replaced by local potentials 
whose strengths are defined in Eq. (||) . In order to preserve this separation of the center-of-mass 
and relative coordinates r and s the interaction must also be separable in r and s. Assuming that 
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the resulting error is small compared to the error introduced by the DME we therefore replace in 
the exchange interaction, as a first approximation, the individual coordinates of the interacting 
particles by their center-off-mass coordinate. This leads to 

v;'{s,v-p)^g^^{v)v:;{s) (53) 

by which we express the exchange interactions ,{s). 



A. Details Of Numerical Calculation 



For our calculations we only consider the stationary ground state of spherical even-even nuclei, 
therefore the single particle wave functions can be separated as 

cj,^ {v) = ?:^^>i,m{e,^) (54) 

r 

The functions '^ijm{0, 0) are spinor spherical harmonics p2| , the index Uq represents the quantum 
numbers = {n, for protons and neutrons, resp. This allows us to simplify the density and 
the kinetic energy density to 



(56) 



The occupation weights for a state aq are determined by pairing correlations. Since it was 
found that pairing gives minor to negligible contributions to the nuclei considered here, the BCS 
approximation was used. A constant pairing matrix element of Gq = 23 MeV/nucleon was used 
and the standard set of BCS equations ||20[] was solved independently for protons and neutrons, 
respectively. 

For protons also the Coulomb interaction contributes to the Hartree-Fock potential U^^ : 

Ucouiir)^ J d^r^pp{r,)j-^^ ~ (^^^ ' e^p^{vfl'' (57) 

The Coulomb exchange term is treated in the local density approximation rather than in the full 
DME because its contribution to the potential energy is small. 

The spin-orbit potential is treated as in ||,^ , leading to a contribution to the potential energy 
Es.o. = J d^r Hs.o.{r) with 

Hs.oXr) = -\wo |p(r) VJ(r) + ^ pq{v) VJ,(r)| . 

The spin-orbit strength Wq is related to the two-body spin-orbit potential in the short-range limit 

by 



Jvls{s)sUs (58) 



Using the M3Y parameterization |2l) one finds a value of Wq ^ 105 MeV fm^ but in correspondence 
to most Skyrme interactions Wq — 120 McV fm^ from reference is chosen. After variation this 
term gives a contribution 



f/-°-(r) = --Wo [p{v) VJ(r) + p,(r) VJ,(r)] , (59) 



to the Hartree-Fock potential and a single particle spin-orbit potential 
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K(^)-lWo^i [p(r)+p,(r)]l.s 



where the spin density is defined as 

-l^T. (2j„, + 1) 



ia,(j«, + 1) - + 1) 



(60) 



(61) 



B. Results for Stable Nuclei 

The various contributions to the nuclear part of the self-consistent DME mean-field in "'"Ca are 
shown in Fig. The isoscalar Hartree potential is weak and repulsive and binding is essentially 
mediated by the Fock parts. The Fock potential is strongly attractive with a depth at the center 
of about -110 MeV. The sum of the Hartree and Fock terms leads to an average potential depth of 
-70 MeV. The direct and exchange rearrangement potentials reduce the depth especially close to 
the center. The mean-field agrees qualitatively with other non-relativistic HF results, e.g [Q. It 
also resembles closely the Schroedinger-equivalent potential from a relativistic density dependent 
Hartree calculation [ p7[p^ ] except for a more pronounced shell structure. 

Of special interest is the influence of the rearrangement potential whose influence has proven 
to be important as found in comparable non-relativistic |43| and relativistic calculations p7| , p^ . 
Rearrangement is of importance for the single particle potential and the energy spectrum. Because 
the rearrangement potential is always repulsive it lowers the single particle potential and leads to 
a change in the single particle spectra and the charge radii. In the calculation of the total binding 
energy the rearrangement effects cancel out. Taking this into account the total binding energy is 
given as 



(62) 



q=p,n 



q—p,n 



with the rearrangement energy calculated from Eq.(27) and Eq.(E8[): 



ni* (r) m 



(63) 
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FIG. 7. Single particle potential for Ca calculated with the D3Y with so = 5. In the lower part, 
contributions to the total isoscalar potential from the direct Uq and the exchange term Uq of the normal 
Hartree-Fock potential as well as the rearrangement potentials (f/g'^, Uq'^) are shown. In the upper 
part, the total isoscalar potential C/q"' as well as its contributions, the attractive HF potential Uo and the 
repulsive rearrangement potential C/(f , are shown. 

An important contribution to mean-field dynamics is provided by the momentum dependence of 
the interaction. Since the exchange parts of the D3Y interaction were found to be most important 
for nuclear binding the question arises whether the momentum structure is described realistically. 
The effective mass m*, Eq.(p7[) is completely determined by the Fock terms, since they are the 
only source of non-localities in the present model. From Fig. || it is seen that a very reasonable 
value of about m*/m ~0.68 is obtained in the interior of a heavy nucleus like ^°*Pb which closely 
agrees with other theoretical results and determinations from proton-nucleus elastic scattering 
. It is worthwhile to emphasize that this result was obtained without readjustments of param- 
eters. The D3Y interaction (or better to say the original M3Y G-matrix |gl|]) has apparently a 
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realistic mixture of short and long range contributions accounting reasonably well for the intrinsic 
momentum structure of an in-medium HF interaction. The results of this section also confirm the 
DME approach. Especially promising results are obtained in conjunction with the effective "Fermi 
momentum" q^, Eq.(^), by which non-local effects are accessible in a self-consistent way. 




r [fm] r [fm] 

FIG. 8. Effective mass of protons (solid line) and neutrons (dashed line) in closed shell nuclei calculated 
with the density dependent D3Y interaction (with the parameterization so ~ 5.) 



TABLE III. Charge radii rc [fm] for closed shell nuclei for different parameterizations of the D3Y 
interaction 





16 Q 


«Ca 


48 Ca 


208pb 


so = 1 


2.79 


3.51 


3.52 


5.51 


So = 5 


2.88 


3.61 


3.60 


5.59 


exp. 


2.73 


3.49 


3.47 


5.50 
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FIG. 9. Charge density distributions of closed shell nuclei from density dependent Hartree-Fock calcu- 
lations using the D3Y interaction. Results for so = 1 (dashed line) and sq = 5 (dotted line) are shown. 
Experimental densities and their uncertainties are denoted by solid lines (taken from Ref. [43]) . 



Fig. |9| compares the calculated charge density distributions of the above closed-shell nuclei with 
experimental data from elastic electron scatterin g p5| . The theoretical point particle distribution 
was folded with a Gaussian proton form factor |43| with {r^)p = 0.8/to. Tab. [II shows the 
calculated charge radii in comparison to the experimental data (taken from p8[|). 

The charge density in the nuclear interior is slightly overestimated for heavy nuclei like ^°*Pb 
and slightly underestimated for the lighter nuclei. In all cases, the density at large radii is higher 
than the experimental data which leads to a overestimation of charge radii as can be seen in Tab. 



III. In general the agreement with the experimental charge radii is satisfying and is, except for 
W), < 2% for the D3Y with sq = 1 and < 4% for the D3Y with sq = 5. In contrast to this, 
the binding energies calculated with sq = 1 are strongly underestimated by more than 0.5 MeV 
per nucleon, whereas the D3Y with so = 5 reproduces the experimental data very well (Tab. 



IV). One sees that an improvement in the binding energies leads to higher values of the charge 



radii. In general it is not possible to reproduce simultaneously both the binding energies and the 
charge radii with one parameter set. This is in agreement with calculations in the relativistic 
density dependent field theory which typically give comparable results But generally the 

agreement with the experimental binding energies of the D3Y with sq = 5 is surprisingly good as 
can be seen in Fig. ^ realizing that the number of used parameters is lower than the one in Gogny 
or Skyrme interactions. 

The observed behavior of the D3Y for different parameterizations shows the importance of the 
density dependence of the interaction. Although both parameterizations have the same properties 
in the range of the saturation density in nuclear matter the results for finite nuclei are quite 
different. This behavior is caused by the differences at low densities. In addition, the rearrangement 
terms enhance the effect of the density dependence. 
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FIG. 10. Binding energy per nucleon for closed shell nuclei from density dependent Hartree-Fock 
calculations using the D3Y interaction. Results for sq = 1 (triangles) and sq = 5 (circles) are shown. 
Experimental binding energies are denoted by open squares (taken from Ref. [^). Binding energies are 
shown for the nuclei ^^O, ""Ca, ''^Ca, ^°Zr, "^Sn and ^"**Pb. Lines are drawn to guide the eye. 



TABLE IV. Binding energy per nucleon Eb/A [MeV] for closed shell nuclei for different parameteriza- 
tions of the D3Y interaction 





16o 


«Ca 


*«Ca 


»«Zr 


i^^Sn 


208pb 


So = 1 


-7.07 


-7.74 


-7.91 


-8.00 


-7.77 


-7.32 


So = 5 


-8.18 


-8.62 


-8.69 


-8.66 


-8.38 


-7.88 


exp. 


-7.98 


-8.55 


-8.67 


-8.70 


-8.36 


-7.87 
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C. The Ground States of Sn Isotopes 



The Sn isotopes are of particular interest for nuclear structure and also astrophysical questions 
because of the closure of the Z=50 proton shell. The known isotopes cover the range from the 
proton dripline at ^'^°Sn |0 to the doubly magic "^^^Sn nucleus which is already /3-unstable. Here, 
we are mainly interested to investigate the isovector properties of the D3Y interaction and, as a 
more general aspect, to test an effective interaction, determined from symmetric nuclear matter 
and stable nuclei, in regions far off stability. No attempt is made to optimize the parameters in 
order to obtain a perfect fit of data. The parameter set with So=5 is used. 

Theoretical and experimental binding energies per particle are compared in Fig. The D3Y 
binding energies are seen to be shifted slightly to the left of the data 1 46| but a fair description of 
the overall dependence on mass and asymmetry is obtained. The strongest binding is obtained for 
"'^^"Sn while experimentally the minimum is found at approximately ^^'^Sn. The binding energies 
of the neutron-rich isotopes are reasonably well reproduced but in view of the deviations on the 
neutron-poor side the agreement might be fortuitous. Towards the proton dripline, i.e. N^Z, 
the binding energies are increasingly underestimated by up to about 0.12 MeV per nucleon close 
to -"^"^Sn. This result is to some extent unexpected because stable nuclei with N'-^Z were well 
described. However, also the binding energy of ^"Zr (Fig. |l0| is underestimated the most and 
apparently the effect is enhanced in going from the N=50 to the Z=50 shell. The deviations 
are in a range typical for calculations with a G-matrix in finite nuclei. Fully phenomenological 
interactions like the Skyrme SLy4 force, for which results are also displayed in Fig. |ll|, and the 
Gogny Dl interaction Ml reproduce the Sn binding energies almost perfectly. 
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FIG. 11. Binding energy per nucleon for the tin isotopes from density dependent Hartree-Fock cal- 
culations using the D3Y interaction with so = 5 (triangles) and the Skyrme SLy4 force (cirles). 
Experimental binding energies are denoted by open squares (taken from Ref. |^^). The lines are drawn to 
guide the eye. 



Variations of the scaling parameter Sq resulted in an overall shift of the binding energy curve 
as it might be expected. Improvements on the neutron-poor side are obtained on the expense 
of a worse description of the neutron-rich isotopes. In particular, the minimum was found to be 
rather insensitive on variations of sq. For So=l the whole curve is shifted up by about 0.5 MeV 
per nucleon with essentially the same shape. These observations lead to the conclusion that the 
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parameters a„ themselves should be readjusted on finite nuclei with constraints to nuclear matter 
rather than retaining the nuclear matter values. Taking the point of view that the nuclear matter 
parameterization accounts for the static part of the interaction, where cluster contributions are 
included effectively by the fitting procedure , such a refinement could be related to non-static 
contributions to the in-medium interaction because of the larger polarizability of the nuclear 
surface. Theoretically, this would amount to go beyond the ladder approximation and to include 
also ring diagrams. Such contributions are obviously contained in e.g. the Skyrmc parameters but 
missing in the D3Y interaction. Over the range of the Sn isotopes the polarization contributions 
can be expected to decrease towards which is known to be a rather stiff nucleus. 
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FIG. 12. Charge density distributions (lower part) and neutron density distributions (upper part) of 
the tin isotopes from density dependent Hartree-Fock calculations using the D3Y interaction. Results for 
So = 5 are shown. 
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Neutron density distributions and charge densities for several Sn isotopes are displayed in Fig. 
[l^ . From ^^°Sn to ^^°Sn the charge distributions are slightly reduced in the interior. As seen 
in Fig. this is accompanied by a mild increase of the charge radius by about 5%. A more 
drastic evolution is found for the neutron densities. Beyond ^■^^Sn an extremely thick neutron skin 
is building up leading to a sudden jump in the neutron rnis radii. The neutron skin thickness is 
more clearly visible in Fig. |lj where the difference of the proton and neutron rms radii is shown. 
In ^'*'^Sn a solid layer of neutron matter is found extending by about 0.48 fm beyond the proton 
matter distribution which is more than a factor of two larger than the neutron skin in ^"^^Pb. The 
same phenomenon is predicted also by the SLy4 interaction except for a somewhat smaller neutron 
excess radius. The saturation of the rms values around A=132 is an indication of the double magic 
nature of ^^^Sn. The increase is directly related to the shell structure in the heavy tin isotopes. 
At A=132 the 2f7/2 shell is filled and pairing does not contribute. At larger masses the neutron 
3p-subshells become populated. Weak binding and the low angular momentum barrier allow a 
large extension of the valence wave functions into the exterior. However, the calculations do not 
lead to a "halo" . 




100 110 120 130 140 

A 

FIG. 13. Root mean square radii of the tin isotopes from density dependent Hartree-Fock calculations 
using the D3Y interaction with sq — 5 for neutrons (squares) and protons (triangles) . The lines are drawn 
to guide the eye. 
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FIG. 14. Difference of the root mean square radii of neutrons and protons of the tin isotopes for the 
D3Y interaction with sq — 5 (triangles) and the Skyrme SLy4 (cirles). The hues are drawn to guide the 
eye. 



An interesting observation is made from Fig. |Ij on the neutron-poor side. Both the D3Y and 
the SLy4 interactions predict a proton skin for A<106. The origins of the proton and neutron 
skins are quite different. Atthe neutron rich side the excess neutrons become less bound because 
of the increasing repulsion from the isovector potential. Since states with i >1 are involved they 
are stored predominantly in the surface region of the nucleus. In the heavy isotopes the isovector 
potential acts attractively in the proton sector and thus balances to a large extent the Coulomb 
repulsion. Since for N^Z the isovector part of the mean-field becomes strongly suppressed the 
repulsion is no longer compensated and Coulomb effects become visible. Hence, the proton skin in 
the neutron-poor Sn nuclides is caused solely by the Coulomb interaction. 

The calculations were extended further into the /3-unstable region. Up to A=160 the neutron 
dripline is still not reached. The separation energies are small but decrease rather slowly. 



V. SUMMARY, DISCUSSION AND CONCLUSIONS 

The density matrix expansion (DME) method was used to derive a local energy density functional 
for non-relativistic meson exchange interactions. The interaction energy density was separated 
into direct and exchange contributions. The Hartree contributions are treated exactly. From the 
systematic expansion of the one-body density matrix provided by the DME an average momentum 
density was determined such that the next to leading order terms in the DME series are canceled. 
In uniform systems like infinite nuclear matter all higher order terms are canceled and the exact 
energy density is recovered. In finite nuclei the approach corresponds to a generalized Slater 
expansion of the density matrix. Hartree-Fock equations were derived. With our choice for q|, 
a self-consistently generated effective mass was obtained which is completely determined by the 
intrinsic momentum structure of the underlying interaction. 

A semi-microscopic approach to the in-medium NN interaction was discussed. The M3Y G- 
matrix from the Paris NN potential was modified by introducing density dependent vertex renor- 
malizations in the isoscalar and isovector particle-hole channels. The functional form of the vertex 
functions was chosen in close analogy to the dynamical structure and the medium dependence 
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of the vertices in a Brueckner G-matrix. The parameters, however, were determined empirically 
by a fit to the binding energy, equilibrium density, compressibility and the symmetry energy of 
symmetric nuclear matter. The momentum structure of the original G-matrix was retained. 

The probably first attempt to parameterize a G-matrix in terms of fixed propagators and density 
dependent vertices goes back to Sprung and Banerjee psf . They approximated the momentum 
structure of the singlet/triplet even-odd interactions obtained from the Reid NN potential by a 
superposition of 5 arbitrarily chosen gaussians with channel dependent in -medi um vertices. The 



functional form of the vertices was chosen in the same way as in section [II A but using always 
a first order polynomial. The nuclear matter binding energy curve of the full Reid G-matrix 
was reproduced the best for a linear dependence of the vertices on kp, corresponding to /3=l/3 
in Eq.(p6|). From the description of individual matrix elements a slight preference for a lower 
value (3—1/6 was deduced However, nuclear matter was found to be badly underbound by 

about 5 MeV since only two-body correlations were taken into account. The renormalized vertices 
obtained in our semi-microscopic approach include effectively also contributions beyond the two- 
body cluster approximation, as was pointed out some time ago by Bethe |p7|l . 

In nuclear matter the D3Y interaction leads to an equation of state with an unusual enhancement 
of attraction at low densities which is not found for phenomenological forces. Unfortunately, G- 
matrix calculations are not available at such low densities but obviously the free NN T-matrix 
should be approached asymptotically. Phenomenological interactions and also the unrenormalized 
M3Y G-matrix do not account for this transition. Without putting too much emphasis on the 
details of this particular result it might be expected that such a region of enhanced attraction 
could indeed exist in the far tails of nuclear wave functions. Before drawing final conclusions on 
this subject a refined treatment of the low density region is clearly necessary. For that purpose an 
extended version of the schematic model for the in-medium vertices might prove to be useful. 

The ground state properties of stable nuclei were reproduced reasonably well. The good descrip- 
tion of binding energies gives confidence to the DME approach, especially in view of the fact that 
an interaction from nuclear matter was used. 
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APPENDIX A: REPRESENTATIONS OF THE NN INTERACTION 

In literature several representations of the NN interaction exist which can be found in, e. g. 
. Here we will only consider the central part of the interaction. The M3Y interaction is given 
in a representation specifying the two-body channel spin and the parity (SE,TE,SO,TO) which is 
commonly used in inelastic scattering. In our calculations we are especially interested in the p-p 
resp. p-n interaction for calculating the Hartree-Fock matrix elements, the (pp, pn) representation. 
It can be easily obtained with the transformation Eof 



V'^ = 
pp 


\ WsE 


+ 61/to] 


nji 






V'^ = 

pn 


\ WsE 


+ Wte^ 


hlVso 


+ SVto] 


np 


pp 


\ WsE 


~ &Vto] 


= V 

nn 






*^ pn 


\ msE 


+ Wte - 


-IVso 


- 3Fto] 


= v 

np 



(Al) 

Here, the indices d resp. e stand for the direct resp. exchange part of the NN interaction. A 
representation especially useful if one uses a contact interaction or is interested in the direct part 
of the interaction is the set (Vq, V^r, Vr, Krr), which follows from the set (SE,TE,SO,TO) through 

Vo = ^ [5VsE + 3Vte + lVso + QVto] 

= 4 i^^^SE + IVte ~ IVso + iVro] 
Id 
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Vr = ^ [IVsE - 3Vte - IVso + SVto] 

V^r = ^[-lVsE-lVTE + lVso + lVTo] (A2) 

Our aim is to find a spin averaged interaction which separates the isoscalar and the isovector 
interaction for both the direct and the exchange part. Defining 

•^0 PP P^ ^P / 

Vr^' = liV^t - - Vnl" + V^t) (A3) 

we find in the (SE,TE,SO,TO) representation 

= ^ [WsE + ^Vte + IVso + Wto\ = Vo 
= ^ [IVsE - Wte - IVso + SVto] = Vr 

V^ = ^ \ZVsE + - IVso - 9yTo] 

V^ = ^ [IVsE - 3Vte + IVso - 3Fto] • (A4) 

One sees that the defined direct terms of this representation are equivalent to the ones of the set 
{Vo, Vcr,Vr,Var) but also that there is no equivalent for the exchange terms. 



[1] W. Schwab, H. Geissel, H. Lenske, K.-H. Behr, A. Briinlc, K. Burkhard, H. Irnich, T. Kobayashi, G. 
Kraus, A. Magel, G. Miinzenberg, F. Nickel, K. Riisager, C. Scheidenberger, B. M. Sherill, T. Suzuki, 
B. Voss, Z.Phys. A350 (1995) 283. 

[2] M. Zinsor, H. Emling, H. Lcnskc ct al., NucI.Phys A (in print); N. Ostrowski. H. G. Bohlcn, A. S. 
Demyanova, B. Gebauer, R. Kapakchieva, Ch. Langner, H. Lenske, M. von Lucke-Petsch, W. von 
Oertzen, A. A. Ogloblin, Y. E. Penionzhkevich, M. Wilpert and Th. Wilpert, Z.Phys. A343 (1992) 
389; H. G. Bohlen, B. Gebauer, M. von Lucke-Petsch, W. von Oertzen, A. N. Ostrowski and H. Lenske, 
Z.Phys. A344 (1993) 381. 

[3] T. H. R. Skyrmo, Nucl. Phys. 9 (1959) 615. 

[4] K. A. Brueckner, J.R. Buchler, S. Jorna and R.J. Lombard, Phys. Rev. 171(1968) 1188. 
[5] H. S. Koehler, Phys. Rep. 18C (1975) 219. 

[6] R. K. Tripathi, A. Faessler, and H. Miither, Phys. Rev. CIO (1974) 2080. 

[7] H. Kiimmel, K. H. Liihrmann and J. G. Zabolitzky, Phys. Rep. 36 (1978) 1. 

[8] J. W. Negelo and D. Vautherin, Phys. Rev. C5 (1972) 1472. 

[9] J. W. Negele, Rev. Mod. Phys. 54 (1982) 913. 
[10] F. T. Baker, L. Bombot, C. DjalaU, C. Glashauser, H. Lenske, W. G. Love, M. Morlet, E. Tomasi- 

Gustafson, J. Van der Wide, J. Wambach, A. Willis, Phys. Rep. (in print). 
[11] M. R. Anastasio, L. S. Cclcnza, W. S. Pong and C. M. Shakin, Phys. Rep. 100 (1983) 327. 
[12] C. J. Horowitz and B. D. Scrot, Nucl. Phys. A464 (1987) 613. 
[13] B. ter Haar and R. Malfliet, Phys. Rep. 149 (1987) 207. 
[14] H. Miither, R. Machleidt and R. Brockmann, Phys. Lett. B202 (1988) 483. 
[15] R. Brockmann and R. Machleidt, Phys. Rev. C42 (1990) 1965. 
[16] R. Brockmann and H. Toki, Phys. Rev. Lett. 68 (1992) 3408. 
[17] H. Lenske and C. Fuchs, Phys. Lett. B345 (1995) 355. 
[18] C. Fuchs, H. Lenske and H. H. Wolter, Phys. Rev. C52 (1995) 3043. 
[19] J. Decharge and D. Gogny, Phys. Rev. C21 (1980) 1568. 

[20] M. A. Preston and R. K. Bhaduri, Structure of the Nucleus, (Addison-Wesley, Inc., 1975). 

[21] N. Anantaraman, H. Toki and G. Bertsch, Nucl. Phys. A398 (1983) 279. 

[22] G. Bertsch, J. Borysowicz, H. McMagnus and W. G. Love, Nucl. Phys. A284 (1977) 399. 



28 



[23] G. R. Satchler, Direct Nuclear Reactions, (Oxford University Press, New York, 1983). 

[24] D. T. Khoa and W. von Oertzen, Phys. Lett. B304 (1993) 8. 

[25] F. J. Eckle, H. Lenske et al., Phys. Rev. C39 1662 (1989); Nucl. Phys. A506 159 (1990). 

[26] H. M. Sommermann, K. F. Ratcliff, T. T. S. Kuo, Nucl.Phys A406 (1983) 109. 

[27] H. A. Bethe, Phys.Rev. 167 (1968) 879. 

[28] H. S. Kochlcr, Phys. Rev. 18C (1975) 217. 

[29] B. D. Day, Rev. Mod. Phys. 39 (1967) 719. 

[30] F. dc Jong and H. Lenske, Phys.Rev. C, in print. 

[31] G. R. Satchler and W. G. Love, Phys. Rep. 55 (1979) 183. 

[32] A. M. Kobos, B. A. Brown, R. Lindsay and G. R. Satchler, Nucl. Phys. A425 (1984) 205. 

[33] S. Haddad and M. Wcigcl, Phys. Rev. C48 (1993) 2740. 

[34] G.R. Satchler, Nucl.Phys. A55 (1964) 1. 

[35] A. L. Fetter and J. D. Walecka, Quantum Theory of Many-Particle Systems, (McGraw-Hill, New York, 
1971). 

[36] E. Chabanat, Doctorial Thesis, University of Lyon (France) 1995; E. Chabanat, P. Bonche, P. Haensel, 

J. Meyer, R. Schaeffer, in print. 

[37] H. Kucharek, P. Ring and P. Schuck, Z. Phys. A334 (1989) 119. 

[38] H. Kucharek, P. Ring, P. Schuck, R. Bengtsson and M. Giro, Phys. Lett. B216 (1989) 249. 

[39] S. Takahara, N. Onishi and N. Tajima, Phys. Lett. B331 (1994) 261. 

[40] J. M. C. Chen, J. W. Clark, E. Krotscheck and R. A. Smith, Nucl. Phys. A451 (1986) 509. 

[41] R. C. Kennedy, Phys. Rev. 144 (1966) 804. 

[42] A.R. Edmonds, Angular Momentum in Quantum Mechanics, (Princeton University Press, Princeton, 
1957). 

[43] D. Vautherin und D. M. Brink, Phys. Rev. C5 (1972) 626. 

[44] C. Mahaux and R. Sartor, Nucl. Phys. A493 (1989) 157. 

[45] P. Grange and M. A. Preston, Nucl. Phys. A219 (1974) 266.; C. J. Horrowitz and B. D. Serot, ihid., 

A368 (1981) 503, A399 (1983) 529 

[46] G. Audi and A.H. Wapstra, Nucl.Phys. A595 (1995) 409. 

[47] M. Chartier ot al., Phys. Rev. Lett. 77 (1996) 2400. 

[48] D. W. L. Sprung, P. K. Banerjee, Nucl.Phys. A168 (1971) 273. 

[49] J. W. Negele, Phys.Rev. CI (1970) 1260. 



29 



